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Abstract 


When using eigenspace assignment to design an aircraft flight control 
system , one must first develop a model of the plant. Certain questions 
arise when creating this model as to which dynamics of the plant- 
need to be included in the model and which dynamics can be left 
out or approximated. The answers to these questions are important 
because a poor choice can lead to closed-loop dynamics that are un- 
predicted by the design model . To alleviate this problem, a method 
has been developed for predicting the effect of not including certain 
dynamics in the design model on the final closed-loop eigenspace. This 
development provides insight as to which characteristics of unmodeled 
dynamics will ultimately affect the closed-loop rigid-body dynamics. 
What results from this insight is a guide for eigcnstructure control law 
designers to aid them in determining which dynamics need or do not 
need to be included and a new way to include these dynamics in the 
flight control system design model to achieve a required accuracy in the 
closed-loop rigid-body dynamics. The method is illustrated for a lateral- 
directional flight control system design using eigenspace assignment for 
the NASA High Alpha Research Vehicle (HARV). 


Introduction 

Fidelity of the design model is a chief concern 
in any control law design process. In this con- 
text, fidelity of the design model corresponds to 
how well a control law, which is designed using this 
model, achieves design objectives when applied to 
the actual system. A fidelity issue was raised in 
the design of a research lateral-directional control 
law (ref. 1) for NASA’s F/A-18 High Alpha Re- 
search Vehicle (HARV). The control law was syn- 
thesized with the CRAFT (control power, robust- 
ness, agility, and flying-qualities trade-offs) design 
methodology (ref. 2), which is a graphical method 
that uses eigenspace placement methods (refs. 3 
and 4). Only rigid-body dynamics were considered; 
other dynamics, such as actuators, were neglected. 
It was noted, however, that when the other dynam- 
ics were included, the closed-loop rigid-pole locations 
varied from those predicted by the low-order design 
model. In this report, a method is developed to de- 
termine which dynamics need to be included in a 
control law synthesis procedure that uses eigenspace 
assignment. 

Elements of the aircraft rigid-body eigenspace are 
well understood, and much is known about desirable 
dynamic characteristics (ref. 5). However, the air- 
craft has dynamics besides those of the rigid body 
(e.g., actuators, control system filters, and transport 
delays). Exactly which dynamics will significantly 
affect the design is normally not known at the outset 
of the control law design process. 


A controller can be synthesized using a system 
with dynamics beyond those of the rigid body, but 
this occurs at a cost. First, the relationship be- 
tween the desired eigenspace and the dynamics of the 
closed-loop aircraft becomes less obvious. Second, 
the speed at which a given set of feedback gains can 
be generated for an iteration of the desired eigenspace 
is reduced. This reduction results in a trade-off be- 
tween the simplicity and speed of a flight control sys- 
tem design iteration and the accuracy of the final de- 
sign. The fundamental question is: what error w r ill 
result if certain dynamics are neglected? The answer 
to this question will provide insight into the relation- 
ship between given unmodelcd dynamics and their 
effect on the rigid-body dynamics of the full-order, 
closed- loop system. 

To provide a clear exposition of the key results, 
the report is organized as follows. The HARV control 
law design is presented in detail as the motivator of 
this research. The rigid-body system used to synthe- 
size the controller and the unmodeled higher order 
dynamics are described. Next, a single- input, single- 
output (SISO) example of the effect of unmodeled dy- 
namics is presented. This example provides the con- 
ceptual basis for the multiple-input, multiple-output 
(MIMO) work presented. 

Symbols 

A stability matrix 

B control matrix 



E 

eigenspace transformation matrix 

a 

angle of attack, deg 

e 

number of effectors 

0 

sideslip angle, rad 

G 

feedback gain matrix 

7 

flight path angle, rad 

9 

gravitational acceleration, g units, 

S 

control effector deflection angle, deg 


1 g — 32.2 ft/sec 2 

C 

damping ratio 

I 

identity matrix 

A 

eigenvalue matrix 

K 

feed-forward gain matrix 

A 

eigenvalue or pole 

k 

reduction factor 

V 

eigenvector 

k F B 

L v ,p,r,b{ j 

feedback gain 

roll moment dimensional stability and 
control derivatives 

T 

4> S 

IV 

time constant 

bank angle, deg 

natural frequency, rad/sec 

l 

number of measurements 

Subscripts: 

M 

measurement matrix 

ail 

aileron 

m 

number of cont rols 

as 

asymmetric stabilator 

N 

feed-through matrix 

b 

body-axis measurement 

N v ,p,r } 8( } 

yaw moment dimensional stability and 

c 

commanded 

control derivatives 

comp 

computed 

n 

number of states of low-order system 

den 

denominator 

n y 

lateral acceleration, g units 

e 

associated with lag in effectors 

P 

roll angular rate, rad/sec 

F 

unmodeled filters in parallel 

Q 

feed-through term effect on augmenta- 
tion matrix 

i 

FB 

unmodeled filter 
feedback 

R 

real number set 

FO 

full-order, closed- loop system 

r 

yaw angular rate, rad /sec 

HO 

unmodeled, higher order dynamics 

s 

T 

TF(s) 

l d 

Laplace transform variable 
time constant matrix 
transfer function 
transport delay time, sec 

i 

l 

LO 

m 

num 

index 

associated with lag in measurement 
low-order, closed-loop system 
associated with lag in control 
numerator 

u 

control vector 

pilot 

pilot commanded 

V 

eigenvector matrix 

rb 

rigid-body pole 

Vt 

airspeed, ft /see 

Orb 

rigid-body dynamics of ( ) 

V 

higher order states 

RO 

roll-off filter 

V s 

lateral velocity, ft/sec 

rtv 

roll thrust vectoring 

w 

unmodcled filter output 

rud 

rudder 

X 

state vector of rigid-body system 

s 

state 

^u,P,<5( j 

side force dimensional stability and 

sens 

sensed value 

control derivatives 

slow 

increased time constant matrix 

z 

2 

measurement vector used for feedback 

ytv 

yaw thrust vectoring 



Superscripts: 

T transpose 

time derivative 
predicted 

* conjugate transpose 

Abbreviations: 

CRAFT control power, robustness, agility, and 
flying-qualities trade-offs 

DEA direct eigenspace assignment 

HARV High Alpha Research Vehicle 

H ATP High- A ngle-of- At tack Technology 

Program 


The low-order, open- loop aircraft model can be 
written as 

x — Ax + Bu 
z = Mx + Nu 

The state, measurement, and control effector vectors 
are 

x T = {v a ,Ps,r a ,<f>8) 

= {pbi r bi n ysensifi comp) 

= (^aib ^rud> ^ytv? ^rtv) 

The measurement equation z is defined by M and 
N matrices to distinguish it from the traditional 
output equation. These measurements are assumed 
to have no noise. The elements of stability, control, 
measurement, and feed-through matrices, at a single 
flight condition, are defined as 


MIMO multiple-input/multiple-output 

SISO single-input / single-output 

Background 

HARV Description 

The analysis presented is motivated by a research 
lateral-directional flight control system design for 
the High Alpha Research Vehicle (HARV), which is 
shown in figure 1. The HARV is part of the NASA 
High-Angle-of- Attack Technology Program (HATP), 
and it will provide flight validation of HATP re- 
search and technology. The HARV is a preproduc- 
tion F/A-18 that has been modified with a thrust 
vectoring system, as shown in figure 2. The thrust 
vectoring is designed to provide additional control 
moments for high angle-of-attack flight. The HARV 
has a research flight control system designed to sim- 
plify the installation and the modification of control 
laws. One intent is to provide flight validation of 
experimental high angle-of-attack control systems. 

Low-Order Aircraft Model 

The research lateral-directional flight control sys- 
tem is designed using linear models of the aircraft at 
various flight conditions. For these models, the rigid- 
body dynamics are fourth order. The states include 
lateral velocity, roll rate, yaw rate, and bank angle. 
The control effectors are aileron, rudder, asymmet- 
ric stabilator, yaw thrust vectoring, and roll thrust 
vectoring. The measurements are roll rate, yaw rate, 
lateral acceleration, and computed sideslip rate. The 
lateral acceleration sensor is located near the pilot 
station, thus preventing the similarity with sensed 
sideslip rate that would occur if it were at the air- 
craft center of gravity. 


" y v yp —vt g cos 7 ' 

L v Lp L r 0 

N v N p N r 0 

0 1 tan 7 0 

^ail ^Aud ^ as \ tv Y S rtv 1 

^^ail ^Aud ^^ytv ^^rtv 

^ail N trud N *as N 6ytv ^rtv 

0 0 0 0 0 _ 

0 cos a — sin a 0 

0 sin a cos a 0 

n Vv n Vp n Vr 0 

8 V 0 p 0 r 0 <j> - 

0 0 0 0 O' 

0 0 0 0 0 

n VS & u n ^ ru(1 n ys as n y*> ytv rtv 

^Vud ^<5ytv ^>rtv - 

Flight conditions range from angles of attack of 
2.5° to 60° at a constant altitude of 25 000 ft in 
unaccelerated flight. 

For this HARV control law design, the states 
chosen above lead to classically defined spiral, roll, 
and Dutch roll modes (at low angle of attack). The 
models used for all examples presented in this work 
are listed in appendix A. 

Low-Order, Closed-Loop System 

The components of the flight cont rol system corre- 
sponding to the low r -order, closed-loop system, shown 
in figure 3, consist of feed-forward and feedback 
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Figure 1. High Alpha Research Vehicle (HARV). 
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Figure 2. HARV thrust-vectoring system. 
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Feed-forward Low-order 

gains plant 



Figure 3. Low-order, closed-loop system. 


gains. The feed-forward gain matrix maps two in- 
puts into the five aircraft control effectors. These 
two inputs are commanded roll and yaw angular ac- 
celerations such that 


ul = ( Pc , i~c) 


The feed-forward gain matrix is a Jacobian of a con- 
trol mapping algorithm, which is discussed in refer- 
ence 6. For the purposes of the research presented, it 
is assumed that the feed-forward gain matrix is given. 
The feedback gains will be used to place closed-loop 
dynamics. Here, the commanded angular accelera- 
tion is referred to as the controls, and the five inputs 
to the aircraft are called effectors. The controls then 
become the sum of pilot input and feedback. The 
feedback gain matrix maps the four measurements 
into these two controls. 

The control system synthesis technique used to 
generate the lateral-directional flight control sys- 
tem was the CRAFT approach based on the di- 
rect eigenspace assignment (DEA) of references 2, 
3, and 4. The CRAFT process provides a graphical 
approach to trade agility, robustness, flying quali- 
ties, and control power. This approach utilizes DEA 
to achieve the desired dynamics selected using the 
CRAFT technique. The DEA generates linear mea- 
surement feedback gains as a function of the design 
model and a desired closed-loop eigenspace. Here, 
the design model contains only the rigid-body air- 
craft dynamics. The design model and control law 
can be expressed as 


x — Ax + Bu 
z = Mx + Nu 
u = Ku c 


U c — U FB + Up}i 0 t 


— Gz 


(System dynamics) (1) 
(System measurements) (2) 
(Feed-forward control) (3) 
(Feed-forward control) (4) 
(Feedback control) (5) 


with n states, m controls, l measurements, and e ef- 
fectors, thus making x G R n } u c € S m , z € R l , and 
u € R e , To derive an expression for the closed-loop 
system, equations (4) and (5) can be substituted into 
equation (3) to get 

u. — KGz + Kupjiot 

Equation (2) is then used for z such that 

u. = KGMx + KGNu + Kiipjj 0 ^ 

Solving for u, 

u = (I - KGN) -1 KGMx + (I - KGN)“ 1 Ku pilot 

By using this equation for u in equation (1), the 
closed-loop system can be stated as 

x = [A + B(I - KGN) _1 KGM]x 
+ B(I-KGN)- 1 Ku pilot 

or in shorthand as 


x = A lo x + B LO u p iiot 

As previously discussed, the feed- forward gains de- 
fine the effector blending used, and the feedback 
gains are synthesized to achieve a desired closed-loop 
eigenspace. For the system described, the feedback 
gain matrix will place l eigenvalues. The DEA will 
also exactly place m elements of the l corresponding 
eigenvectors, or alternatively, it will achieve a least- 
squares fit of i elements of an individual eigenvector, 
where m < i < n. 

For the purposes of feedback gain synthesis, there 
are four states, four measurements, five effectors, and 
two controls. This implies exact placement of the 
poles for the low-order, closed-loop system. Also, 
two eigenvector elements can be exactly placed for 
each mode. The feedback gains are designed to be 
scheduled with angle of attack. The feed-forward 
gains are scheduled primarily with angle of attack, 
dynamic pressure, and thrust. Appendix B contains 
baseline feed-forward and feedback gain matrices. 

Full-Order, Closed-Loop System 

The complete aircraft plus control system has 
dynamics beyond that of the system previously de- 
scribed. This includes other aircraft dynamics, such 
as actuator dynamics, and additional elements of 
the control system, such as structural notch filters. 
The layout of this control system and plant model is 
shown in figure 4. 
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Feed -forward Low-order 



T+ Second-order Actuator 


Filters models 

Feedback Second-order 


gams 

filters 

LO 
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Figure 1. Layout of HARV lateral-directional control system 
and plant model which differentiates low-order and higher 
order, unmodeled dynamics. 


Models of the actuator dynamics are available in 
reference 7. These actuator models range from first 
order to eighth order for each of the five actuators. 

Various flight control system filters are also not 
included in the design process. These include first- 
order, roll-off filters of 25 rad/sec, which are placed 
on each of the controls as part of the flight control 
system design. A notch filter is also placed on each of 
the controls as well as on each of the measurements. 
An exception is the sensed lateral acceleration chan- 
nel that has two notch filters in series. 

Here, the HARV full-order, closed-loop system 
will have 25 states. In the following development, 
however, not all unmodeled dynamics will be initially 
considered. The term full-order, closed-loop system 
will be used to describe a closed-loop system con- 
sisting of the open-loop, rigid-body dynamics, the 
feed- forward and feedback gains, and the particular 
unmodeled dynamics under consideration. 


Spectral Decomposition 

To evaluate full-order, closed-loop system dynam- 
ics, consider the portion of the eigenspace that cor- 
responds to the rigid-body system. When the un- 
modeled dynamics are of sufficiently higher frequency 
than the low-order, closed-loop system, or they are 
outside the bandwidth of the pilot, the dynamics that 
are dominant to the pilot will be determined by the 
rigid-body eigenspace. The rigid-body eigenspace 
is a combination of the rigid-body eigenvalues and 
eigenvector elements that are associated with the 
rigid-body modes. As an example, the rigid-body 
roll pole of the low-order, closed-loop system and of 
the 25-state system is shown in figure 5 as a function 
of angle of attack. 


Low-order system roll mode 



a, deg 


Figure 5. Low-order system and full-order, closed- loop system 
roll pole. 


The full-order, closed-loop system has dynamics 
defined by 


x = A F qx + B^QUpjjot 


The rigid-body eigenspace is characterized here by a 
spectral decomposition of A F q such that 

A f0 V = VA 


The eigenvalue and eigenvector matrices are then 
partitioned to separate the eigenspace of interest so 
that 

v= rv„ v 12 - 

.V21 V22. 

A = [ Al °1 

.0 A*2 . 

where Ai is a square matrix containing achieved 
rigid-body eigenvalues on its diagonal and Vji is a 
square matrix containing achieved rigid-body eigen- 
vector elements associated with rigid-body states. 
Ideally, one would like to compare the eigenspace of 
the low-order, closed-loop system with the eigenspace 
of 

( A Fo)r6 = VnAiV^ 1 

(i.e., the rigid-body eigenspace of the full-order, 
closed- loop system). This matrix has the same di- 
mensions as the low-order, closed-loop system matrix 
A LO- 
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Although this reduced-order system is not a sub- 
stitute for a look at frequency response and time 
histories of the full-order system, it does allow the 
control law designer to compare the dynamics of the 
full-order, closed-loop system with well- understood 
aircraft rigid-body dynamics. Also, knowledge of the 
achieved rigid-body eigenspace is important because 
its relationship with the desired dynamics fundamen- 
tally determines whether or not certain dynamics can 
be omitted from the design model. 


Actuator Low-order 

model plant 



Eigenspace Transformation Matrix 

A concept that will appear repeatedly in this work 
is that of a single matrix transformation in the form 

- E ( A Fo)r6 

where A[q is the low-order, closed-loop system ma- 
trix and (A fo)rb * s composed of the rigid-body por- 
tion of the full-order, closed-loop system eigenspace, 
as described in the section entitled “Spectral De- 
composition.” The term E will be referred to 
as an eigenspace transformation matrix. Although 
A iq and (A Fo)rb are n °t actually the eigenspaces, 
their eigenvalues and eigenvectors are the rigid-body 
eigenspace of the low- and full-order closed-loop 
systems, respectively. All three of these matrices 
are square, and they have dimensions equal to the 
number of states in the rigid-body system. As 
the unmodeled dynamics become less significant to 
the achieved rigid-body eigenspace, the eigenspace 
transformation matrix should approach the identity 
matrix. 

A look at the most significant contributors to 
the eigenspace transformation matrix may provide 
insight into what causes the rigid-body dynamics to 
differ from those of the lbw-order, closed-loop system. 
If the effect on the rigid-body eigenspace can be 
easily calculated, such a calculation could be useful 
in determining whether certain dynamics can be left 
out of the design model. 

SISO Example 


Figure 6. SISO example with first-order roll mode approxi- 
mation at angle of attack of 5°. 

For the purposes of control law design, the first- 
order lag aileron actuator is not included in the 
design model. This SISO low-order, closed-loop 
design model can be written as 

M = (6) 

w pilot( 5 ) 5 ^ail^FB 

With the actuator model, this full-order, closed-loop 
system has two states, and it can be written as 

p(«) = r ^i (7) 

u piIot(s) ( T ail s + l)(s - L p ) - L 6tU k FB 


In state space, this full-order, closed-loop system is 
expressed as 



The low-order, closed-loop system has one state, with 
the eigenvalue 


^ LO = L p + L Sgi] k FB (8) 


which is the desired roll mode pole. The full-order, 
closed-loop stability matrix is 


a fo 


L P L tai\ 
T -J k FD ~tF} 


To facilitate the discussion of the multiple-input/ 
multiple-output (MIMO) results in this report, a 
single-input/single-output (SISO) example is now 
presented. The SISO system (fig. 6) presented is a 
first-order roll mode approximation of the HARV at 
an angle of attack of 5°. Aileron deflection is the 
plant input, and roll rate is the plant output. The 
first-order lag aileron actuator model is the unmod- 
eled dynamics in this example. 


Of interest are the eigenvalues and eigenvectors of 
this closed-loop system corresponding to the rigid- 
body dynamics, or roll mode. The eigenvalue is 
A po rb an d the eigenvector is v\\ when this full-order, 
closed-loop system matrix is written as the spectral 
decomposition 




i'll ^12 


i'll ^12 


Q 

cr 

0 

- T J k FB ~rj. 


^V2\ ^22 _ 


.^21 ^22. 


o 

A//o. 
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Note that A fo b and v\\ are scalars and that, in the 
absence of v<i\ , is arbitrary. In this development, 
only the upper left and lower left partitions of the 
above matrix equation are needed such that 

Lpv n + L^v 21 = ^ 11 A FOrb (9) 

r^kpBVn ~ = "2 1 A FO rb (10) 

Now, the goal is to eliminate V 2 \ and solve for 
Aj io i n terms of A FO rb - The method that proves 
successful when working with the MIMO system is 
the following. Equations (9) and (10) are rewritten 
as 

LpV 11 - vu\FO rb = ~ L 6^2\ (11) 

kpBVU - "21 + TuMl^FOrb (12) 

Equation (11) is expanded into an infinite series 
such that each term contains the right-hand side of 
equation (12). To start this expansion into an infinite 
series, add and subtract — L p T a [\iy2i^F(y rb and group 
terms to give 

LpV 11 - "l 1 AfOrfc = - ^,("21 + TgSlV2lXFO r b) 

+ ^-'6 lli | 7 ’aU i/ 21 A/O r (, 

Continuing in a similar manner by adding and 
subtracting 

(~ l Y L S M T ln l/ 2l ApOri, (i - 2,3, , oo) 

and grouping terms to get the right-hand side of 
equation (12) in each term yields the infinite scries 

L p v u - I'll ^FO rb = + Tni!"21 A/.o rb ) 

+ 1-6^1 T a il( "21 + Tail "2 1 A TO ril ) A FO rb 
~ L *ail T ail^ y21 + T ai) "21 A FO rb ) ^FO rb 
+ Z '6 ai i r ail("21 + T Bi l"21 XFO rb )>FO rb • • • 

Then eliminate "21 by using equation (12) so that 
£p" 11 - "llAFO r6 = -Lb^k FB Vu + L ban T a \]kFB>‘FO rb l/ U 

~ L S ai] T l\ k FB>' 2 FO rb l/ U 

+ L/> m\ T ail k FB *FO rb "1 > 

Dividing out the now arbitrary v\\, the above 
becomes 

Lp - A FO rb = -Ls aii k FB + Tail k FB^FO rb 

L b &a T M k FB>?FO rb + ks^k FB Xl 0rb . . . 


One can now relate the low-order (cq. (8)) and full- 
order, closed-loop rigid-body pole as 

A w = X FO rl , + ^ ail T ai ifc Ffl A FOjfj 

“ L 6^ T l\\k F B>^FO rb + L 6 aii T a\) k FB>?FO rb • ’ ’ 
This is an infinite series that converges as long as 

l r ail A FO rb \ < 1 

In the case where the actuator has a much 
lower time constant than the rigid-body mode 
I r ail A FO r & I ^ I? the first two terms constitute a good 
approximation to the infinite series such that 

A LO « (l + ^ a ji T ai 1 ^FB ) A FO rb 

or 

A FO rb ~ A FO rb = ( L + L b lM T ai\kFfj) ^LO 

where A FO rb denotes an approximation of A FO rb ' A 
plot of this approximation is shown in figure 7. Note 
that A FO rb approximates A Fo rb u P an actuator 
time constant of approximately 0.06 sec, which is 
three times slower than the actual actuator time 
constant. 



Figure 7. Effect of first-order actuator model on SISO system 

and approximation. 

Because of the divergence of A 2q and A^J^, a 
point will exist where the open-loop actuator time 
constant is high enough that it becomes necessary 
to include this actuator model in the design process. 
This point will be determined by the required level 
of accuracy in the final design. A point will also 
exist where the rigid-body roll mode and the actuator 
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mode couple together to form a single oscillatory This full-order, closed-loop system now has n + e 
mode. states with dynamics expressed as 


Approximating the Transformation 
First-Order Actuators 

In this section, an approximation of the eigen- 
space transformation is derived and tested for un- 
modeled dynamics placed at the input of the rigid- 
body aircraft. The eigenspace transformation relates 
the rigid-body eigenspace of the low-order, closed- 
loop system used for the design with the full-order, 
closed-loop system. Here, the SISO results shown in 
the section entitled “SISO Example” are generalized 
for first-order lag actuator dynamics applied to the 
MIMO design. 

Without feed-through term . Consider the 

following system, shown in figure 8, which represents 
the HARV rigid-body system and the first-order lag 
actuator models: 

x = Ax + Bv (System dynamics) 
z — Mx (System measurements) 

T € v = — v + u (Actuator dynamics) 

u = KGz + Kupflot (Control law) 




A 

T^KGM 




t- l k 


u 


pilot 



The closed-loop stability matrix for the low-order 
system is 

A L0 = A + BKGM 

and the closed-loop stability matrix for this full-order 
system is 


A FO 


A B ' 

T e -1 KGM -T" 1 . 


Of interest are the eigenvectors and eigenvalues of 
this closed-loop system corresponding to the rigid- 
body dynamics. These correspond to Vp and Ai 
when this full-order system matrix is rewritten as 
the spectral decomposition, 


A B 


< 

< 

to 


< 

< 

to 


> 

O 

i 

T.r'KGM -T ” 1 


_V 2 1 V *22 


V 21 v 22 . 


. 0 A 2 . 


(14) 


Here, we need only the upper left and lower left 
partitions of the above matrix equation, or 

AV u +BV 2 i = VnAi (15) 

T~ 1 KGMVn - T“ 1 V 2 i = V- 2 i A 1 (16) 

or equivalently, 

VnAj - AV n =BV 2 i (17) 


Figure 8. System block diagram with first-order lag actuators 
as unmodeled dynamics. 


The measurement feed-through term is omitted from 
this derivation for clarity. The actuator matrix is 
diagonal and contains the time constant associated 
with each actuator: 


n 


T e = 


0 


0 

T e 


This formulation corresponds to actuator dynamics 
with the transfer functions 


Vj{s) _ 1 

Ui(s) TjS + 1 


(*= l,2,...,c) 


KGMVn - V 2 i + T c V 21 Ai (18) 

As an intermediate step to obtain the eigenspace 
transformation approximation, equation (17) is ex- 
panded as an infinite series. To generate this scries, 
add and subtract BT t V 2 iAi and group terms cor- 
responding to the right-hand side of equation (18). 
This becomes 

V U A! - AVn = B(V 2 i + T c V 2 i AO - BT c V 2 i A x 

Then add and subtract BTjVoi A'f and group as 
before such that 

VnAi - AVn = B(V 21 + T e V 21 A0 - BT e (V 2l 
+ T c V 21 A0Ai + bt2v 21 a? 
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where A 2 = AA, A 3 = AAA, . . . Continuing to 
add and subtract BT*V 2 iA* 1 (i = 3,4,... , oo) and 
following the same grouping strategy yields the fol- 
lowing infinite series: 


The matrix 2-norm of A equals maxA(AA*) 1 / 2 . By 
virtue of being eigenvectors, Vn can be multiplied 
by an arbitrary scalar. This factor is chosen such 
that the matrix norm of KGMVn is unity, so that 


VnAi - AV n = B(V 2i + T C V 21 AO 

- BT e (V 21 + T e V 21 A 1 )A 1 
+ BT2(V 21 +T c V 21 A 1 )Af 
-BT3(V 21 +T, : V2iA 1 )A? + ... 

Equation (18) can now be substituted into each term 
on the right to get 

VnAi - AVn = BKGMVn - BT e KGMV n Ai 

+ BTgKGMVnAi - . . . (19) 

Postmultiply equation (19) by Vn . The infinite 
series that results relates the rigid-body eigenspace 
of the low-order, closed-loop system with that of the 
full-order, closed-loop system in which 

A LO ~ (A Fo)rb + BTeKGM(A/ro) rb 

- BT2KGM(A FO ) r 2 6 

+ BTgKGM(Afo)^ - . . . (20) 


||T*KGMVn A, || < (||T e || ||Ai||)*‘ 

Convergence of equation (19) is therefore guaranteed 
when 

l|T e || || A! || <1 

Feed-through term . For many systems, the low- 
order model of interest may have a feed-through or 
N matrix term. This is the case for the HARV control 
law design. Here, the previous analysis and the 
development are extended to obtain the eigenspace 
transformation matrix for this type system. To 
begin, consider the system that represents an aircraft 
with first-order lag actuator models in which 

x = Ax + Bv (System dynamics) 

z = Mx + Nv (System measurements) 
T e v = —v + u (Actuator dynamics) 

u = KGz + Kupji ot (Control law) 

The system has n + e states, with the closed-loop 
dynamics 


where 

(Afo)rt = VnAiVfj 1 

as defined in the section entitled “Spectral 
Decomposition.” 

As the amount of augmentation increases (K 
and G), the effect of the unmodeled dynamics on the 
rigid-body dynamics is increased. As the actuators 
get faster (he., as the elements of T e approach 0) the 
effect on the closed-loop eigenspace is decreased. The 
eigenspace of the low-order, closed-loop system and 
the rigid-body eigenspace of the full-order, closed- 
loop system become equal. 



A 

T-'KGM t^ 1 


+ 


0 

T-’K 


u pilot 


B 

KGN-T-’ 



z = [M N]j*J 


The spectral decomposition of the system matrix is 


A B 

'Vn 

V 12 

T“'KGM T^KGN-T" 1 . 

.V 2 i 

V 22 . 


The convergence of equation (19) is equivalent to 


V u 

V 12 ‘ 

'Ai 0 ' 

the convergence of 


v 21 

v 22 . 

.0 A 2 . 


^(-lynKGMVnAi 

?=0 

This series can be bounded using the matrix 2-norm 
and the Schwarz inequality (ref. 8) such that 

||T’KGMV n Ai|| < ||T e |P||KGMVn|| HAiH 4 

(i = 0, 1, — , oo) 


where A\ and Vn correspond to the rigid-body 
eigenspace. As before, only the upper left and lower 
left partitions of this equation are used, or 

VnAj-AVn = BV 21 (21) 

KGMV n - Q t V 21 + T e V 2 iAi (22) 

where 

= I - KGN 
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This is a useful definition for the development to 
follow, since from the section entitled “Low-Order, 
Closed-Loop System” 

A l0 = A + BQ“ 1 KGM 

As in the previous section, equation (21) is ex- 
panded as an infinite series with terms grouped to 
match the right-hand side of equation (22). The 
series is generated as before: add and subtract 

BQ c , T e V 21 A 1 to equation (21), so that 

VnAj - AVu = BQ-^QeVai + TW^iAO 
-BQ- 1 T e V 2 iA 1 

Then continue in a similar manner as before and 
apply equation (22) as follows: 

VnAj - AVu = BKGMVn - BQ.'TcQJ’KGMVh Aj 

+ BQ~ 1 T < .Q“ 1 T,.Q” , KGMViiAi - . . . 



Figure 9. 2-norm of terms in infinite series. 


This expression leads to the concept of a single 
matrix (an eigenspace transformation matrix) which 
transforms the low-order, closed-loop stability matrix 
to approximate the rigid-body eigenspace of the full- 
order, closed-loop system such that 


Convergence is guaranteed when [|Q e 1 T e || |[Ai|| < 1. 
To obtain the final form, postmultiply by V^ 1 , so 
that 

Alo = (A^o)rfc + BQ; 1 T r QJ 1 KGM(A^o)r6 

- BQ" 1 T e Q" 1 T C Q“ 1 KGM( A/. o)j! fe + . . . (23) 

An important simplification of equation (23) is 
suggested by multiplying through by (A Fo)rb ' so 
that 

A^AfoU 1 -I + BQ- T T e Q fc - l KGM 

- BQ- 1 T,Q t 7 1 T,Q- 1 KGM(A FO ) r 6 + - . - 

The 2- norm of each successive term on the right-hand 
side, except the identity term, is shown in figure 9. 
Clearly, the first two terms constitute a reasonable 
approximation to the infinite series for this example. 


( A FC>)r6 ~ E lA LO ( 24 ) 

For this full-order, closed-loop system, the eigenspace 
transformation matrix is simply 

E = I + BQ^TcQ-'KGM 

Application 

The results of the previous development are ap- 
plied to the HARV lateral- directional flight control 
system design. The validity of the previous analysis 
in predicting the effect of first-order actuator models 
on the rigid- body dynamics is shown. The low-order 
system wa^used to generate the feedback gains. Be- 
cause there are four measurements, the rigid-body 
poles are exactly placed. Moreover, the low-order, 
closed-loop eigenspace is the desired eigenspace. 

In this example, the actuators are modeled as the 
following transfer functions: 


An advantage of using the first two terms as an 
approximation to the series is the following predic- 
tion of the rigid-body eigenspace of the full-order, 
closed-loop system expressed in terms of the low- 
order, closed-loop system: 

(A FO ) r6 ~ (I + BQ^TeQ^KGMJ^Ato 


r TFa.i\i s ) ~ i 1 

4 S + 1 


TF as (s) — , 

30 s + 1 

TF T t v (s) = -j — 
® s + 1 


^•^rudC 5 ) ~ j . 

30 s + 1 

TFytvis) = , 

35 s + 1 
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The approximate eigenspacc transformation matrix 
is then 



E = I + BQ“ 1 


0 


1 

w 

1 

30 


0 


1 

15 


Q^KGM 


15 J 


(25) 


The expression E _i A/^Q is postulated to be an ap- 
proximation to the rigid-body eigenspace of the full- 
order, closed-loop system (i.e., VnAiV^ 1 ). 


The first check on the validity of the approxima- 
tion will be on the convergence of the infinite series 
used to derive the eigenspace transformation matrix. 
The test 

HQc 'Tell ||Ai < 1 


is applied using the 2-norm. At an angle of attack of 
20°, the left side of the equation is equal to 0.0737. 
This value indicates convergence of the series. 


To examine the range in which the eigenspace 
transformation matrix is valid, the poles of the ac- 
tuators are gradually slowed down. This procedure 
is done by multiplying the actuator poles by a reduc- 
tion factor k which is varied from 1 to 0 such that 


where 


T 


-1 

slow 


kT 


-l 

e 


T s iowV = -v + u 


■o 

c3 


-1.6 i- 
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T3 

a 


c 

3 
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Q 


3.0 r- 




1.0 h 


i 
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The error associated with using this eigenspace trans- 
formation matrix in predicting the roll mode pole 
and Dutch roll natural frequency for a flight condi- 
tion with an angle of attack of 20° is shown in fig- 
ure 10. The actuators have been slowed by a factor 
of 5 when the aileron actuator and the roll mode poles 
form an oscillatory pair. The unmodeled dynamics 
will need to be much closer in frequency to the low- 
order, closed- loop system for there to be a problem 
with the assumptions made. 

Figure 11 shows how well E -1 A iq predicts the 
rigid-body roll mode pole. The low^-order characteris- 
tics are those of the design model, and they represent 
placed dynamics. The full-order characteristics are 
those of the low-order plant plus the actuator dy- 
namics. The predicted characteristics are obtained 
from equation (24), and they arc expected to predict 
the full-order characteristics. The Dutch roll mode 
pole is shown in figure 12. Note that the variation of 
the roll mode and Dutch roll frequencies are captured 
along with the variation in Dutch roll damping. The 


Figure 10. Plots that verify accuracy of eigenspace transfor- 
mation matrix approximation with angle of attack of 20° 



Figure 11. Roll mode pole prediction with actuator models. 
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a, deg 

Figure 12. Dutch roll mode pole prediction with actuator 
models. 


negligible variation of the spiral mode pole, which 
is not shown, is also captured by the approximation. 
For all cases, the approximate eigenspace transforma- 
tion matrix accurately predicts the rigid-body poles 
of the full-order, closed-loop system. 

The predicted effect of actuators on the desired 
eigenvector elements is shown in figure 13. The first 
of these elements is the 0-to-/? ratio in the Dutch roll 
mode. The second element is the /3-to-p ratio in the 
roll mode. The magnitude of the ratio of each of these 
two eigenvector elements is shown. These ratios give 
an indication of how much the roll and Dutch roll 
modes are coupled. The low-order /?-to -p ratio in the 
roll mode is 0 for all cases. The low-order (p-to-ft ratio 
in the Dutch roll mode coincides with that of the 
full-order, closed-loop system. The /3-to-p ratio is in 
seconds, and the (p-to-ft ratio is nondimensional. The 
eigenvectors appear to be affected considerably less 



Figure 13. Eigenvector prediction with actuator models. 


than the eigenvalues by unmodeled dynamics. This 
result will show up repeatedly in the examples that 
are presented. Again, the approximate eigenspace 
transformation matrix accurately predicts the full- 
order, closed-loop rigid-body eigenspace. 

Dynamics in Other Locations 

General Form of Eigenspace 

Transformation Matrix 

One would like a general expression for the 
eigenspace transformation matrix for dynamics at 
other points, including the input to the plant, in the 
control system structure. It is possible to construct 
the system with a first-order lag at three locations 
within this system. The eigenspace transformation 
matrix for a first-order lag in any or all of these loca- 
tions will then be readily available from the general 
form. 
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Feed-forward 


Low-order 



Figure 14. System block diagram with first-order lags in three 
locations. 


Consider the system shown in figure 14, in which 


The eigenvector matrix has been partitioned such 
that the first-column partition contains n columns, 
and the second contains e + m-bZ columns. The first-, 
second-, third-, and fourth-row partitions contain n, 
e, m, and Z rows, respectively. Equations in the first- 
column partition of the matrix equation (26) are 

VnAi-AVn =BV 21 (27) 

KV 31 = V 2 i + T p V 2 i A| (28) 

GV 41 =V 3 i+T ro V 31 A 1 (29) 

MVn = -NV 2 i + V 41 + T / V 41 A 1 (30) 

To get the feed-through term in a useful form, a 
second version of equations (28), (29), and (30) is 
needed, so that 


x = Ax + Bv e T c v c = -v e + u 


KV 31 = Q«V 2 i + T e V 21 A 1 + KGNV 21 (31) 


u — Kv m T m v m — v m + u c 

U FB = Gv I T tvi = -v/ + z 

z = Mx + Nv e u c = u FB + u pi i ot 

where x E R n , v e E R e , v m € R m > and u/ € 
R? . When the equations governing this system are 
combined, the full-order, closed-loop system becomes 


' X ' 


' A 

B 

0 

0 ' 


r X 1 

V r 


0 


T-’K 

0 


V e 


► = 





< 


V m 


0 

0 

_T _i 

A m 

T m 'G 


v,„ 

< Vi > 


_T“ ! M 

Tf l N 

0 

-Tr 1 J 


» v^ > 




0 

0 

T -1 

A in 
0 


/ Upilot 


with n + e + m + l states. The spectral decomposition 
of the system matrix becomes 


' A 

B 

0 

0 ' 


'V u 

Via' 

0 

_ T -i 

e 

T^K 

0 


V21 

V 22 

0 

0 

_ T” 1 

x m 

T^G 


V31 

v 32 

-Tf'M 

T/ _1 N 

0 

-TfL 


.V41 

v 42 . 


■v„ 

V21 

V31 

.V 41 


V12I 


V 22 

V 32 

V42. 


'A! 

0 


0 ' 
a 2 . 


(26) 


where Vy and A, are the eigenvector elements and 
eigenvalues associated with the rigid-body eigenspace 
of the full-order, closed-loop system, respectively. 


G V 41 - Q m V 31 + T to V 31 Aj + GNKV 31 (32) 

MVn = Q/V 4 i + T/V 4 i Ai + NKGV 41 - NV 2! 

(33) 

where 

Q e = I - KGN 
Qm = I - GNK 
Qi = I - NKG 

As in the previous development, expansion ob- 
tained via adding and subtracting needed terms will 
turn equations (27) through (30) into an infinite se- 
ries. It is again assumed that the speed of the un- 
modeled dynamics is faster than that of the low- 
order, closed-loop system. Hence, each new term of 
the series will be smaller than its predecessor. To 
form an analogous approximation to the scries, all 
terms containing two or more time constant matrices 
and two or more eigenvalue matrices will be dropped. 
Expand equation (27) so that equation (31) can be 
used; thus, 


V u Ai - AV„ = BQ e - ! (Q,V 21 + T e V a iA|) 

- BQ“ 1 T,,Q J 1 (Q e V 21 + TeVjiAOAi 
+ BQ” 1 T ( .Q < 7 i T c Q“ i (Q(.V 2 1 
+ T e V 2 iAi)A?-.„ 

Then apply equation (31) and drop the terms with 
higher powers of T e and Aj, so that 


VnA, - AV U « BQ ( -'(KV 31 - KGNV 21 ) 

- BQ e , T c Q7 1 (KV 31 - KGNV 2 i)Ai 


For the remainder of the development, this equation 
will be stated as an approximation; all terms with 
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more than a single appearance of T and will be 
dropped as they occur. 

Expand the above equation again so that the 
right-hand side of equation (29) appears as 

VuA| - AV„ » BQ^K^i +T m V 3 iAi - GNV 2I ) 
-BQ- , KT m V 3 iA l 

- BQ t - , T t Q ( 7 1 K(V 3 i +T m V 3 iAi 
-GNV 2 i)Ai 

Then by equation (29), 

VuAj - AV n « BQ t 7 1 K(GV 4 i - GNV 2l ) 

- BQ''KT„|V 3 iAi 

- BQ ( r 1 T ( Q- , K(GV 4 i - GNV 21 )A] 

Equation (32) will now be used to eliminate the last 
V 31 term, so that 

V n A, - AV u « BQ e ”'KG(V 4 , - NV 21 ) 

- BQ e -’KT m Q- , (Q m V 31 + T m V 31 A 1 )A 1 

- BQ e r'T e Q~ 1 KG(V 4 i - NV 2 i)Ai 

and then 

V„Ai - AV„ « BQ,7 i KG(V4i - NVji) 

- BQ- , KT m Q- , G(V 4 i - NKV 31 )A] 
-BQ^TVQ^KGCVTn -NVaiJAi 

Equation (28) is then used; thus, 

V„Ai - AV n « BQ e -*KG(V4i - NV 21 ) 

- BQ t 7 ] KT m Q- 1 G(V4i - NVjiJAi 

- BQ“ 1 T e Q" , KG(V 4 i - NV 21 )Ai 

Expand the above equation again so that the 
right-hand side, of equation (30) will allow V 41 to 
be eliminated: 

VuAi-AVn « BQ“ 1 KG(V 4i + T/V 41 Ai 

-NV 21 )-BQ- 1 KGTiV 4 iAi 

- BQ;’KT m Q m 1 G(V 4 i + T;V 4 i A] 

- NV«)Ai - BQ<r 1 TeQ~ 1 KG(V 4 j 
+ T/V 41 Ai - NV 2 i)Ai 


Then apply equation (30), so that 

VnA, - AV 11 » BQ-^GMVn 

— BQ“ 1 KGT;V 4 iAi 

- BQJ 1 KT w Q“ 1 GMV 11 Ai 

- BQ^TeQ-'KGMVnAi 

Next expand the above expression so that equa- 
tion (33) is used and 

Vn A, - AV U ~ BQj’KGMVn 

-BQ- 1 KGT,Q,- , (V 4 i + T]V4iA 1 )A 1 
-BQ c 1 KT m Q- , GMV„Ai 
- BQ-’T. Q-’KGMVn Ai 

or 

VnA] - AV n « BQ^KGMV,, 

- BQ‘ 1 KGT ( Qj" 1 (MVn - NKGV 41 

+ NV 2 i)A] -BQ f - , KT m Q- 1 GMV u A, 

- BQ“’T P Q ( 7 , KGMV 11 Ai 

Equations (29) and (28) are subsequently applied to 
eliminate the last V 4 j term; thus, 

VnA] - AVn ~ BQ“ ] KGMV U 

- BQ, 'KGTjQ-'fMVn - NV 2 , 

+ NV 2 i)A] - BQ“ 1 KT m Q“ 1 GMVn A] 

- BQ" 1 T ( Q t 7 1 KGMV] 1 Ai 

All unwanted terms are eliminated, so that 

Vn A, AV n « BQ^KGMVn 

- BQ J 1 KGT; Q ; 1 MV 1 ] A] 

- BQ ( T 1 KT m Q” 1 GMVn A] 
-BQj’TeQ-’KGMVnA] 

Finally, note that 

Q-’KG = (I - KGN)-'KG = KG (I NKG )- 1 = KGQ ' 1 
Q“’K = (I - KGN)“'K = K(I GNK ) -1 = KQ ” 1 
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by the matrix inversion lemma. The approximation 
is then converted to the familiar form: 


VnAj - AVu « BQ^KGMVn 

- BKGQ^T/Qj ^MVnA] 

- BKQ- 1 T m Q^ 1 GMV n A 1 

- BQ-^eQ-'KGMVnAi 

Postinultiplying the result by yields 

A-lo ~ (Afo)r6 + BKGQ ; 'T;Q ( ] M(ApQ) r ft 
+ BKQ m T m Q w GM(A/?o) r j 

+ BQ-^eQ-^GMfA^oU 

where 

A w = A + BQ-’KGM 

(A FO ) r6 = V 11 A 1 V 1 - 1 1 


Application 

By using the subset of equation (34) that corre- 
sponds to a system with first-order lag in controls, 
one obtains 


E = I + BKQ m 1 T m Q“ 1 GM (35) 


where 

Qm = I - GNK 

Alo = a + bq-’kgm 

These results will be used to predict the effect of the 
roll-off filters shown in figure 4. The roll-off filter 
transfer function is 


TF ro {s) 


1 

+ 1 


in each channel. Therefore, for this example, the 
eigcnspace transformation matrix is 


E = I + bkq - 1 


l 

25 


1 

25 


Qm GM 


(36) 


The eigenspace transformation matrix is now 

E = I + BKGQ^TjQ-’M + BKQ “ 1 TmQ” 1 GM 
+ BQ~ 1 T e Q" 1 KGM (34) 

Note that the previous derivation for first-order 
lag actuator models is a subset of this result. If 
the first-order lag at any one of these points is not 
present, then the associated time constant is 0. In 
addition, within each matrix of time constants, any 
of the elements can be 0 when that particular signal 
does not have a lag to be modeled. 

Two or more first-order lags in series at a partic- 
ular location in the system are handled in a simple 
way. The time constants associated with each chan- 
nel only need to be summed before applying equa- 
tion (34). For example, the approximate eigenspace 
transformation matrix of two first-order lags in series 
with time constants of 0.01 sec and 0.02 sec 

( : 1 V 1 ) = I 

Vo.Ols + 1 / \0.02s + 1 J 0.002s 2 + 0.03s + 1 

1 

~ 0.03s + 1 

will be the same as that of a single first-order lag 
with a time constant of 0.03 sec. Note that this 
substitution will be good over the frequency range 
of the rigid-body dynamics. 


The purpose of the roll-off filter is to attenuate 
high-frequency commands. This attenuation should 
reduce potential problems caused by structural dy- 
namics or noise. Results similar to those of the actu- 
ators are expected, since both the actuators and the 
roll-off filters are first-order lags. Figure 15 shows 
the predicted effect of the filters on the roll mode 
pole. The predicted effect of the filters on the Dutch 
roll mode pole is shown in figure 16. Figure 17 
shows the predicted effect of the filters on modal 
coupling. The effect of these unmodeled dynamics 
is more pronounced because the roll-off filters are 
slower than the actuators. Again, the approximate 
eigenspace transformation matrix accurately predicts 
the rigid-body eigenspace of this full-order, closed- 
loop system. 

Higher-Order Filters 

Equivalent First-Order Lag 

Now that the eigenspace transformation is avail- 
able to predict the effects that first-order lags have 
on the rigid-body eigenspace of the full-order, closed- 
loop system, it is useful to generalize the transfor- 
mation to higher order filters. More specifically, one 
would like to handle dynamics such as lead or lag 
compensators, transport delays, and notch filters. 
Handling these dynamics is essential when faced with 
the second-order filters found in the HARV flight con- 
trol system. A system with unmodcled higher order 
filters on the controls is shown in figure 18. The feed- 
through term has been removed for this analysis. 
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Figure 15. Roll mode pole prediction with roll-off filters. 




Figure 17. Eigenvector prediction with roll-off filters. 


Command 



Figure 16. Dutch roll mode pole prediction with roll-off 
filters. 


Figure 18. System block diagram with second-order command 
filters as unmodeled dynamics. 
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The system can be written as follows: 

x = Ax + BKw (System dynamics/ 

feed-forward gains) 

z = Mx (System measurements) 

u F £ = Gz (Feedback control) 

u c = u FB + Upiiot (Feed- forward control) 

v = A F v + B F u c (Filter system dynamics) 

w = M F v + N pu c (Filter measurements) 

where the filters in parallel channels are realized as 



thus giving each higher order filter a transfer function 
of the form 

My. (si - AyJ' 1 By. +Ny (* = 1,2, . . . ,m) 

The closed-loop system matrix spectral decompo- 
sition for this system is 

A + BKN f GM BKM f 1 r Vi i V 12 " 
B F GM A f J [ V 2 i v 22 | 

= “Vn V 12 l I" Ai 0 - 

.V21 V22. . 0 a 2 . 

where again Ai and Vp correspond to the rigid- 
body eigenspace. When multiplied out, this equation 
yields four matrix equations. Two of these equations, 
which correspond to the upper left and lower left 
partitions, can be written as 

VnAi-(A+BKN F GM)Vn = BKM f V 2 i (37) 


-Af’BfGMVn = V 21 - A^VaiAi (38) 

where it is assumed that the higher order filters have 
no poles at the origin. Equation (37) is expanded by 
methods used in previous sections to yield the infinite 
series 

V n A! - AV n -BKN F GMVn 

= BKM F (V 2 1 - A~W 2 iA x ) + BKM F A“ 1 (V 2 1 
-A- 1 V 2 iA,)A 1 +BKM;.'A- 2 (V 21 -A- , v 21 A 1 )Ap ... 

Substituting equation (38) yields 

Vn A| - AV n - BKNFGMVn 
= -BKM F A~ 1 B F GMV il 

- BKM F A“ 2 B F GMV n Aj - BKM F A“ 3 GMV U A 2 X - ... 

Postmultiplying by V^ 1 , we get 

(A F o)r6 = A + BKN f GM - BKM F A F 1 B F GM 

- BKM F A^ 2 B f GM(A f o) r6 

- BKM f A f 3 BfGM(A TO ) 2 6 - . . . (39) 

The type of higher order filters considered is 
restricted to those with a steady-state gain of 1 
through all channels. The steady-state gain is the 
limit of the filter transfer function as s goes to 0. 
This gain of 1 yields the relations 

- My,A^ lB y. +Ny. = 1 

or 

— M f A^*B f T N f = I 

Using this result to combine the first three terms on 
the right-hand side of equation (39), we get 

Alo = {^Fo)rb + BKM F A i7 2 B F GM(A F (o) r ^ 

+ BKM f A- 3 B f GM(A fo )^ + ... (40) 

The approximate eigenspace transformation matrix 
of any higher order filter with a steady-state gain of 
1 is the same as that of a first-order lag with an equiv- 
alent time constant matrix equal to M F A^ 2 B F . 

Second- order filter, A second-order filter mod- 
eled as the transfer function 

Wi(s) _ ^den-j / S 2 + 2(num^num,S + ^num ; \ 
'Ui(s) \ “h ^Cden^denj 5 F CJ 2 eii . J 
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has the following system matrices 


or in state-space form 


A/, 

0 — ^den; 


1,2,... 

,m) 

B /, 

deiij 

r / i 2 

^num, 

~ ^den,; 



i “ u;2 

^ruirrij 

2 (Cmim^niimj 

i — Cden, 

^deiij) _ 


M fi 

= [0 1] 




N /, 

den, 




u/ num i 





The equivalent time constant is MyvAy 2 By. or 


A /,= 


-i 


^cierLj- 


M/i — 1 


B h = 


T den ? ; ~ r numj 
r denj 


N/, = 


^rnmij 

^deiij 


The rigid-body eigenspace of the full-order, closed- 
loop system that includes this filter can be accurately 
predicted using the eigenspace transformation ma- 
trix corresponding to a first-order lag with a time 
constant of 

T i = r den* — T minij 

This result implies that such a compensator could be 
designed to cancel much of the effect of unmodeled 
dynamics on the rigid-body eigenspace. 


n = 2 


Cde 


k ^den, 


w num ; 


(41) 


Application 

The HARV lateral-directional second-order filters 
have the transfer function 


As an example, consider a notch filter with a de- 
nominator damping ratio of 0.7, a numerator damp- 
ing ratio of 0.1, and equal natural frequencies in 
the numerator and the denominator. From equa- 
tion (41), we see that the effect of the filter on the 
rigid-body eigenspace will be virtually the same as 
that of a first-order lag with a break frequency of 
83 percent of the natural frequency of the notch filter. 

Time delay . The first-order Pade approxima- 
tion to a transport delay is written as 


s 2 T 2£ mim cj num s + ^num j 
s 2 + 2ClerMienS + ^den ) 

with characteristics shown in table 1. The first two 
characteristics listed are the command filters; the re- 
maining five are measurement notch filters. The two 
lateral acceleration filters are configured in series. 
Each of these two filter sets will now be considered. 

Table 1. HARV Second-Order Filters 


TF(s) 


u 2 


Wj{s) _ 1 - (<d,/ 2 ) 5 

u ii s ) l+(^d,/2)s 
The system matrices can be written as 


A /,= 


-2 

t d i 



M fi = 1 N f. = — 1 


Filtered 

signals 

^den i 

rad/sec 

Cdcn 

^niim y 

rad /sec 

Cnum 

Pc 

40 

0.6 

140 

0.74 

r c 

40 

.6 

140 

.74 

Pb 

80 

.7 

80 

.08 

rb 

150 

.7 

150 

.08 

hysons 

58 

.7 

58 

.08 


80 

.7 

80 

.08 

Pscns 

80 

.7 

80 

.10 


The appropriate first-order lag for this transfer func- 
tion is found to have a time constant equal to the 
transport delay magnitude such that 


For the command filters, equations (34) and (41) 
imply that the eigenspace transformation matrix 
should be 


T i — tdi 


E = I + BKQ m ] 


0.0194 

0 


0 

0.0194 


Q^GM (42) 


Lead/lag compensator . Another type of filter 
that has been considered is a lead/lag compensator 
of the form 

mis) _ T nuni; ,s + i 

u *( s ) ^en^ 1 


Although a feed-through term was not included when 
deriving the equivalent time constants for higher 
order filters, it has been postulated that the feed- 
through term should be handled in the same way as 
suggested in previous sections. 
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Dutch roll damping ratio Dutch roll natural frequency, rad/sec 



Figure 19. Roll mode pole prediction with command second- 
order filters. 




Figure 20. Dutch roll mode pole prediction with command 
second-order filters. 



Figure 21. Eigenvector prediction with command second- 
order filters. 


From the previous discussion, these second-order 
command filters, having zeros with frequencies of 
140 rad/sec and poles with frequencies of 40 rad/sec, 
will affect the rigid-body cigcnspace of the full-order, 
closed-loop system in the same manner as a roll- 
off filter with a time constant of 0.0194 see. The 
eigenspace results are shown in figures 19, 20, and 21. 
Again, the approximate eigenspace transformation 
matrix accurately predicts the resulting rigid-body 
eigenspace. 

The second-order filters in the measurement loop 
are notch filters designed to cancel resonant peaks 
in the structural model. There are two second- 
order filters in series for the n yscns channel, so the 
approximate time constant associated with each will 
be summed to approximate the effect of both. The 
eigenspace transformation matrix associated with the 


20 



Low order 

With measurement notch 



Figure 22. Roll mode pole prediction with measurement notch 
filters. 



Figure 23. Dutch roll mode pole prediction with measurement 
notch filters. 



Figure 24. Eigenvector prediction with measurement notch 
filters. 


measurement notch filters is, from equations (34) 
and (41), 


E = I+BKGQ, 1 


'0.0155 

0 


0.0083 

0.0369 


0 


Qf’M 


0.0150 J 

(43) 


The effects of these notch filters are shown in fig- 
ures 22, 23, and 24. The eigenvalues and eigenvectors 
ofE -1 A LO predict very well the resulting rigid-body 
eigenspace. Only a slight degradation in the predic- 
tion of the small /3-to-p ratio in the roll mode effect 
is observed. Hence, the approximate first-order lags 
are a suitable substitute for the higher order filters 
in this example. 

This section has shown how the effects of higher 
order filters on the closed-loop rigid-body eigenspace 
can be evaluated using first-order lags placed at the 
same loop location. This result is expected as long as 


21 



the steady-state gains of the unmodeled higher order 
filters are unity. 

HARV Full-Order, Closed-Loop System 

An approach is now proposed to predict the si- 
multaneous effect of all unmodeled dynamics (which 
is introduced in the section entitled “Full-Order, 
Closed-Loop System”) on the HARV control law de- 
sign. The strategy is straightforward: approximate 
the HARV full-order, closed-loop system of figure 4 
using the closed-loop system in figure 14 by replac- 
ing all unmodeled dynamics with equivalent first- 
order lag components. Equation (34), which defines 
the eigenspacc transformation matrix in this simpler 
case, is then directly applied. 



The full-order, closed-loop system is placed in the 
simpler form as follows. Note that the actuators are 
already modeled by first-order lags; therefore, noth- 
ing needs to be done. Turning to the other un- 
modeled dynamics, the second-order command and 
measurement notch filters are replaced by approxi- 
mate first-order lags derived in the section entitled 
“Higher Order Filters.” Here, it is postulated that 
the results of that section, which were only consid- 
ered filters at a single-loop location, extend to filters 
at multiple-loop locations. Next, the filters in series 
(i.e., the first-order roll-off and approximate com- 
mand filters) are replaced by single first-order lags. 
The time constant of the approximate lag equals the 
sum of the time constants of the replaced lags; this 
leads to a system of the form shown in figure 14 
with first-order lag elements at three separate loop 
locations. Equation (34) can then be used to obtain 
the eigenspace transformation matrix, where (from 
eq. (25)), 


'0.0208 


0.0250 


T e = 


0.0333 


0.0208 


0.0208 


From equations (36) and (42), 


'0.04 

0 


'0.0194 

0 


'0.0594 

0 



+ 



~ 



0 

0.04 


0 

0.0194 _ 


0 

0.0594 


From equation (43), 


T, 


0.0155 

0 


0.0083 

0.0369 


0 ' 
0.0150. 


Figure 25. Roll inode pole prediction with full-order model. 



Figure 25 shows the prediction of the roll mode 
pole. Figure 26 shows the Dutch roll mode pole. 


Figure 26. Dutch roll mode pole prediction with full-order 
model. 
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Figure 27. Eigenvector prediction with full-order model. 


Figure 27 contains the eigenvector element ratio pre- 
diction. As before, the eigenvalues and the eigen- 
vectors ofE^A^o predict the rigid-body eigenspace 
of the 25th-order, closed-loop system. 

Conclusions and Recommendations 

The approach presented in this report allows the 
control law designer to predict to what extent un- 
modeled dynamics affect the closed-loop, rigid-body 
eigenspace. Such insight is important when decid- 
ing which dynamics will be included in the model 
used to design the controller. A single-input, single- 
output example was used to illustrate and predict the 
effect of unmodeled dynamics on a closed-loop, rigid- 
body pole. This result was extended to multiple- 
input, multiple-output dynamics, thus leading to 
the concept of an eigenspace transformation ma- 
trix that relates the desired closed-loop rigid-body 


eigenspace to that obtained in the presence of un- 
modeled dynamics. 

The approach was first developed for unmodeled 
first-order lag elements at one specific loop location 
and then extended to multiple-loop locations. The 
eigenspace transformation matrix was shown to ac- 
curately predict the achieved rigid-body eigenspace 
for this type of unmodeled dynamics. For higher 
order unmodeled filters with a steady-state gain of 
unity, derived approximate first-order lag compo- 
nents were shown to be a suitable replacement in 
predicting the achieved rigid-body eigenspace. The 
approximate components were easily found from the 
state-space expressions of the unmodeled higher or- 
der filters. Also, the aggregate effect of many types 
of unmodeled dynamics on the achieved rigid-body 
eigenspace was shown to be well predicted. 

In conclusion, note that there will always be 
some errors in the achieved rigid-body eigenspace, 
whether they are caused by poor mathematical mod- 
els, off-design flight conditions, or unmodeled dynam- 
ics. The goal is to find what these errors are and 
determine if they will be acceptable. A method for 
predicting how some unmodeled dynamics affect the 
rigid-body eigenspace is now available. 

Some important areas of future work suggested 
by this research include the following: 

1. Developing a method to predict the achieved 
rigid-body dynamics when a significant frequency 
separation does not exist between the rigid-body 
and unmodeled dynamics (i.e., when the convergence 
criteria would be violated) 

2. Formulating sensitivity relationships that de- 
termine changes in the achieved rigid-body dynamics 
caused by incremental changes in unmodeled equiv- 
alent first-order lag time constants 

3. Developing a desired eigenspace adjustment 
procedure to account for the effect of unmodeled dy- 
namics on the final full-order, closed-loop eigenspace 
using the eigenspace transformation matrix concept 

4. Developing methods to change the design 
model, instead of the desired eigenspace, to account 
for the effect of unmodeled dynamics in the control 
system design process 

NASA Langley Research Center 
Hampton, VA 23681-0001 
January 31, 1994 
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Appendix A 

HARV Lateral-Directional Aircraft 
Model 

This appendix contains the linear models repre- 
senting the High Alpha Research Vehicle (HARV) 
lateral-directional rigid-body aircraft dynamics at 
the 13 flight conditions used in this study. A full 
nonlinear model (ref. 7), written in the Advanced 
Control Simulation Language (ACSL), was used to 
generate these Jacobians. The form of the linear sys- 
tem is 

x = Ax + Bu (Low-order dynamics) 

z = Mx + Nu (System measurements) 

with four states , four measurements , and five effec- 
tors. The states are lateral velocity, stability-axis roll 
rate, stability-axis yaw rate, and bank angle (given in 
units of feet per second, radians per second, radians 
per second, and radians, respectively). The measure- 
ments are body-axis roll rate, body-axis yaw rate, lat- 
eral acceleration, and sideslip rate (given in units of 
radians per second, radians per second, g units, and 
radians per second, respectively). The effectors are 
aileron, rudder, asymmetric stabilator, yaw thrust 
vectoring, and roll thrust vectoring (all given in units 
of degrees). 

Finite differencing was used to generate the 
Jacobians. Table Al contains the perturbation step 
size on states and effectors used to generate the 
Jacobians. These Jacobians are listed in table A2. 


Table Al. Finite Differencing Perturbation Sizes 


State 

Perturbation 

size 

Effector 

Perturbation 
size, deg 

v s 

Vr/lO ft/sec 

^ail 



Vs 

.08 rad/sec 

Aiid 



r s 

.08 rad/sec 




4>s 

.04 rad 

^ytv 





^rtv 

j 



The open-loop roll and spiral mode poles for each 
flight condition are shown in figure Al. Figure A2 
shows the open-loop Dutch roll frequency and damp- 
ing. Figure A3 contains plots of the two open-loop 
eigenvector characteristics used in the numerical ex- 
amples to quantify modal coupling. Note that the 
0- to-/3 ratio is nondimensional. 



Figure Al. Roll and spiral mode poles of open-loop aircraft. 


Natural frequency, rad/sec 



<j>-to*P ratio in Dutch roll 


B-to -p ratio in roll mode, sec 



Figure A3. 4>- to-/3 in Dutch roll and (3-to-p roll mode eigen- 
vector element ratios of open-loop aircraft. 


24 




Table A2. HARV Jacobi ans 


a = 2.5° 

V T = 837 ft/sec 




A - 

-0.1855 

0.0072 

-836.1471 

32.1675 


-0.0271 

-2.7580 

0.6086 

0.0000 


0.0058 

0.1084 

-0.1614 

0.0000 


0.0000 

1.0000 

0.0000 

0.0000 


B = 

-0.0946 

0.5304 

-0.2000 

0.2628 

0.0056 

0.4398 

0.0578 

0.4538 

-0.0048 

0.0181 

-0.0209 

-0.0419 

-0.0171 

-0.0317 

-0.0015 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.9990 

-0.0436 

0.0000 


0.0000 

0.0436 

0.9990 

0.0000 


-0.0030 

0.0964 

-0.0560 

0.0000 


-0.0002 

0.0000 

-0.9993 

0.0384 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0197 

-0.0009 

-0.0218 

-0.0041 

-0.0008 

-0.0001 
a = 5° 

0.0006 

Vp — 598 ft /sec 

-0.0002 

0.0003 

0.0000 

A - 

-0.1305 

0.1512 

-597.5921 

32.1675 


-0.0187 

-1.5272 

0.6757 

0.0000 


0.0050 

0.1152 

-0.1529 

0.0000 


0.0000 

1.0000 

0.0000 

0.0000 


B = 

-0.0551 

0.2975 

-0.1025 

0.2461 

0.0024 

0.2746 

0.0314 

0.2224 

-0.0059 

0.0176 

-0.0283 

. -0.0242 

-0.0161 

-0.0294 

-0.0019 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.9962 

-0.0872 

0.0000 


0.0000 

0.0872 

0.9962 

0.0000 


-0.0021 

0.0535 

-0.0462 

0.0000 


-0.0002 

0.0003 

-0.9992 

0.0538 
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Table A2. Continued 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0134 

0.0002 

-0.0100 

-0.0038 

-0.0007 

-0.0001 
a = 10° 

0.0005 

Vt ~ 421 ft/sec 

-0.0002 

0.0004 

0.0000 

A = 

-0.0955 

0.1610 

-420.2861 

32.1674 


-0.0194 

-0.8687 

0.6852 

0.0000 


0.0062 

0.1333 

-0.1871 

0.0000 


0.0000 

1.0000 

0.0000 

0.0000 


B - 

0.0286 

0.1409 

-0.0358 

0.2500 

0.0006 

0.1350 

0.0146 

0.1010 

-0.0085 

0.0177 

-0.0261 

-0.0130 

-0.0161 

-0.0291 

-0.0032 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.9848 

-0.1736 

0.0000 


0.0000 

0.1736 

0.9848 

0.0000 


-0.0012 

0.0295 

-0.0392 

0.0000 


-0.0002 

0.0004 

-0.9991 

0.0765 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0068 

-0.0002 

0.0042 

-0.0038 

-0.0007 

-0.0001 
a = 15° 

0.0003 

V T = 361 ft/sec 

-0.0001 

0.0006 

0.0000 

A = 

-0.0702 

0.1331 

-360.9948 

32.1679 


-0.0185 

-0.5227 

0.6580 

0.0000 


0.0069 

0.1218 

-0.2358 

0.0000 


0.0000 

1.0000 

0.0000 

0.0000 


B = 

-0.0198 

0.0873 

-0.0149 

0.2497 

-0.0013 

0.0817 

0.0078 

0.0618 

-0.0105 

0.0174 

-0.0236 

-0.0089 

-0.0160 

-0.0282 

-0.0045 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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Table A2. Continued 


M = 

0.0000 

0.9659 

-0.2588 

0.0000 


0.0000 

0.2588 

0.9659 

0.0000 


-0.0007 

0.0169 

-0.0407 

0.0000 


-0.0002 

0.0004 

-0.9994 

0.0891 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0044 

-0.0002 

-0.0026 

-0.0037 

-0.0006 

-0.0001 
a — 20° 

0.0002 

Vt = 334 ft /sec 

0.0000 

0.0007 

0.0000 

A = 

-0.0558 

0.0399 

-333.9384 

32.1655 


—0.0236 

-0.2994 

0.6024 

0.0000 


0.0095 

0.0932 

-0.2796 

0.0000 


0.0000 

1.0000 

0.0000 

0.0000 


B = 

-0.0076 

0.0559 

-0.0079 

0.2406 

0.0038 

0.0463 

0.0037 

0.0451 

-0.0121 

0.0160 

-0.0183 

-0.0061 

-0.0168 

-0.0262 

-0.0063 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.9397 

-0.3420 

0.0000 


0.0000 

0.3420 

0.9397 

0.0000 


-0.0005 

0.0069 

-0.0477 

0.0000 


-0.0002 

0.0001 

-1.0002 

0.0963 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0026 

-0.0002 

-0.0021 

-0.0036 

-0.0007 

0.0000 
a = 25° 

0.0002 

Vp = 307 ft/sec 

0.0000 

0.0007 

0.0000 

A = 

-0.0472 

-0.0569 

-308.1857 

32.1649 


-0.0191 

-0.2336 

0.5856 

0.0000 


0.0103 

0.0901 

-0.3405 

0.0000 


0.0000 

1.0000 

0.0000 

0.0000 
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Table A 2. Continued 


B - 

0.0108 

0.0375 

0.0061 

0.2548 

0.0027 

0.0315 

0.0005 

0.0355 

-0.0158 

0.0169 

-0.0162 

-0.0035 

-0.0180 

-0.0260 

-0.0084 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.9063 

-0.4226 

0.0000 


0.0000 

0.4226 

0.9063 

0.0000 


-0.0002 

0.0007 

-0.0541 

0.0000 


-0.0002 

0.0002 

-1.0006 

0.1044 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0015 

-0.0001 

-0.0018 

-0.0037 

-0.0008 

0.0000 
a = 30° 

0.0001 

Vy = 282 ft /sec 

0.0000 

0.0008 

0.0000 

A - 

-0.0403 

-0.1336 

-282.2590 

32.1549 


-0.0099 

-0.3858 

0.7811 

0.0000 


0.0060 

0.2001 

-0.5262 

0.0000 


0.0000 

1.0000 

-0.0280 

0.0000 


B - 

0.0195 

0.0270 

0.0183 

0.2471 

0.0007 

0.0211 

-0.0006 

0.0288 

-0.0167 

0.0166 

-0.0137 

-0.0022 

-0.0185 

-0.0243 

-0.0099 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.8660 

-0.5000 

0.0000 


0.0000 

0.5000 

0.8660 

0.0000 


-0.0007 

0.0041 

-0.0675 

0.0000 


-0.0001 

-0.0005 

-1.0009 

0.1140 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0008 

0.0000 

-0.0013 

-0.0036 

-0.0008 

0.0001 

0.0001 

0.0001 

0.0009 

0.0000 
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Table A2. Continued 


a = 35 ° Vp = 268 ft /sec 


A = 

- 0.0423 

- 0.2074 

- 267.5458 

31.9377 


- 0.0027 

- 0.3022 

0.8106 

0.0000 


0.0019 

0.1766 

- 0.6487 

0.0000 


0.0000 

1.0000 

0.1202 

0.0000 


B = 

0.0246 

0.0214 

0.0250 

0.2432 

0.0007 

0.0155 

- 0.0011 

0.0246 

- 0.0183 

0.0155 

- 0.0124 

- 0.0015 

- 0.0196 

- 0.0226 

- 0.0112 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.8192 

- 0.5736 

0.0000 


0.0000 

0.5736 

0.8192 

0.0000 


- 0.0012 

- 0.0048 

- 0.0715 

0.0000 


- 0.0002 

- 0.0008 

- 1.0010 

0.1195 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

- 0.0005 

- 0.0001 

- 0.0011 

- 0.0036 

- 0.0008 

0.0001 
cv = 40 ° 

0.0001 

Vp = 261 ft/sec 

0.0001 

0.0009 

0.0000 

A = 

- 0.0435 

- 0.2705 

- 261.6818 

31.4210 


0.0003 

- 0.3069 

0.6522 

0.0000 


- 0.0018 

0.2127 

- 0.6251 

0.0000 


0.0000 

1.0000 

- 0.2193 

0.0000 


B = 

0.0257 

0.0181 

0.0277 

0.2401 

0.0007 

0.0125 

- 0.0018 

0.0210 

- 0.0198 

0.0143 

- 0.0123 

- 0.0008 

- 0.0208 

- 0.0208 

- 0.0124 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.7660 

- 0.6428 

0.0000 


0.0000 

0.6428 

0.7660 

0.0000 


- 0.0018 

- 0.0081 

- 0.0628 

0.0000 


- 0.0002 

- 0.0010 

- 1.0008 

0.1202 
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Table A2. Continued 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0004 

-0.0001 

-0.0011 

-0.0036 

-0.0008 

0.0001 
a = 45° 

0.0001 

Vr = 262 ft/sec 

0.0001 

0.0009 

0.0000 

A - 

-0.0383 

-0.2061 

-261.3775 

30.5801 


-0.0105 

0.1843 

0.0830 

0.0000 


0.0088 

-0.2481 

-0.1543 

0.0000 


0.0000 

1.0000 

-0.3264 

0.0000 


B = 

0.0079 

0.0162 

0.0265 

0.2349 

0.0007 

0.0109 

-0.0030 

0.0247 

-0.0211 

0.0131 

-0.0129 

0.0004 

-0.0269 

-0.0189 

-0.0135 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.7071 

-0.7071 

0.0000 


0.0000 

0.7071 

0.7071 

0.0000 


-0.0012 

0.0350 

-0.0210 

0.0000 


-0.0001 

-0.0008 

-0.9994 

0.1169 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0009 

-0.0001 

-0.0011 

-0.0036 

-0.0008 

0.0000 
a = 50° 

0.0001 

V T = 262 ft/sec 

0.0001 

0.0009 

0.0000 

A = 

-0.0333 

0.0412 

-261.3253 

29.4804 


-0.0077 

-0.0330 

0.2272 

0.0000 


0.0084 

-0.0162 

-0.3423 

0.0000 


0.0000 

1.0000 

-0.4366 

0.0000 


B = 

-0.0118 

0.0147 

0.0367 

0.2294 

0.0007 

0.0087 

-0.0035 

0.0240 

-0.0222 

0.0117 

-0.0129 

0.0014 

-0.0290 

-0.0170 

-0.0144 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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Table A 2. Continued 


M = 

0.0000 

0.6428 

-0.7660 

0.0000 


0.0000 

0.7660 

0.6428 

0.0000 


-0.0008 

-0.0122 

-0.0208 

0.0000 


-0.0001 

0.0002 

-0.9985 

0.1126 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0016 

-0.0001 

-0.0003 

-0.0036 

-0.0008 

0.0000 
Q = 55° 

0.0001 

Vj = 266 ft/sec 

0.0001 

0.0009 

0.0000 

A - 

-0.0380 

0.1644 

-265.4073 

28.4186 


-0.0065 

-0.0508 

0.2138 

0.0000 


0.0087 

0.0197 

-0.3784 

0.0000 


0.0000 

1.0000 

-0.5303 

0.0000 


B = 

-0.0145 

0.0143 

0.0719 

0.2260 

0.0007 

0.0064 

-0.0030 

0.0309 

-0.0229 

0.0103 

-0.0124 

0.0011 

-0.0380 

-0.0150 

-0.0152 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M = 

0.0000 

0.5736 

-0.8192 

0.0000 


0.0000 

0.8192 

0.5736 

0.0000 


-0.0009 

-0.0050 

-0.0240 

0.0000 


-0.0001 

0.0006 

-0.9990 

0.1070 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0017 

0.0002 

0.0018 

-0.0036 

-0.0008 

-0.0001 
a = 60° 

0.0001 

Vt = 276 ft/sec 

0.0003 

0.0009 

0.0000 

A = 

-0.0472 

0.1523 

-275.5784 

26.6075 


-0.0058 

-0.0638 

0.2364 

0.0000 


0.0089 

0.0320 

-0.3957 

0.0000 


0.0000 

1.0000 

-0.6794 

0.0000 
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Table A2. Concluded 


B = 

—0.0156 

0.0144 

0.1012 

0.2251 

-0.0071 

0.0051 

-0.0022 

0.0227 

-0.0236 

0.0107 

-0.0131 

-0.0001 

-0.0296 

-0.0128 

-0.0144 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

M 

0.0000 

0.5000 

-0.8660 

0.0000 


0.0000 

0.8GG0 

0.5000 

0.0000 


-0.0013 

-0.0083 

-0.0097 

0.0000 


-0.0002 

0.0006 

-0.9995 

0.0965 


N = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0018 

-0.0003 

0.0037 

-0.0034 

-0.0001 

-0.0001 

0.0001 

0.0004 

0.0008 

0.0000 
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Appendix B 

Baseline HARV Flight Control System 
Gains 

This appendix contains the linear feed-forward 
and measurement feedback gains for the 13 flight 
conditions considered. There are two controls and 
five effectors. As previously stated, the effectors are 
aileron, rudder, asymmetric stabilator, yaw thrust 
vectoring, and roll thrust vectoring (which arc all 
measured in degrees). The measurements are body- 
axis roll and yaw rate (which are measured in radians 
per second), sensed lateral acceleration (which is 
measured in g units), and sensed sideslip rate (which 
is measured in radians per second). The controls 
are commanded roll and yaw acceleration (which are 
measured in radians per second squared). The feed- 
forward and measurement feedback gain matrices are 
listed in table Bl. 

The feed- forward gains are designed such that the 
two controls are mapped into the five effectors. The 


two controls are commanded roll and yaw angular 
acceleration. The feed- forward gain matrix is a 
Jacobian of a control mapping algorithm (discussed 
in ref. 6), which is evaluated at the various trimmed 
flight conditions. For the purposes of this research, 
the set of feed-forward gain matrices is assumed to 
be given. The controls are a sum of pilot input and 
feedback (fig. 3), 

u = K(u fb + u p j] ot ) 

The measurement feedback gain matrices map the 
four measurements into the two controls such that 

U FB ~ Gz 

The baseline gains listed w r crc designed using the 
CRAFT (control power, robustness, agility, and 
flying-qualities trade-offs) procedure of reference 2. 
These gains are designed to be scheduled with angle 
of attack. 
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Table Bl. Feed- Forward and Feedback Gain Matrices 


a = 2.5° 


K = 




2.4595 

1.6767 



-1.5926 

-34.1549 



0.3513 

2.5021 



0.0000 

0.0000 



0.0000 

0.0000 



G = 




-0.1062 

-0.6541 

0.0441 

-0.2824 

-0.0642 

0.1223 

0.0460 

2.1799 

a = 5° 




K = 




4.0206 

2.4478 



-5.6765 

-62.0159 



0.8725 

4.9350 



0.0000 

0.0000 



0.0000 

0.0000 



G = 




-0.5442 

-0.7925 

-0.0019 

-0.3825 

-0.0625 

0.1133 

0.0524 

1.7372 

a = 10° 




K = 




6.9486 

2.5159 



-21.1127 

-116.4603 



2.3793 

8.9879 



0.0000 

0.0000 



0.0000 

0.0000 



G = 




-1.1081 

-1.2024 

-0.1014 

-0.1852 

-0.0437 

0.1778 

0.0607 

1.7539 


a = 15° 

K - 

10.1457 -0.9134 

-27.6594 -103.2218 

3.4755 7.2567 

-6.1403 -21.9256 

0.0000 0.0000 
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Table BL Continued 


G - 


-1.2663 

-1.1728 

-0.0005 

-0.6436 

0.0228 

0.2459 

0.0600 

1.5931 

P 

II 

to 

o 




K = 




14.2474 

-5.3125 



-33.5753 

-90.8832 



5.0302 

3.8308 



11.3867 

-30.2217 



0.0000 

0.0000 



G = 




-1.4859 

-2.0383 

-0.6063 

-1.2320 

0.0681 

0.3996 

0.2004 

1.9451 

a = 25° 




K = 




14.1097 

-14.4480 



-41.1839 

-84.1229 



4.7607 

-1.7777 



-13.5142 

-28.1494 



0.0000 

0.0000 



G = 




-1.5959 

-1.9635 

0.3622 

-1.1257 

0.0885 

0.4970 

0.0183 

2.1152 

a = 30° 




K = 




13.2511 

-23.2949 



-47.7514 

-75.2633 



3.9250 

-7.5768 



-15.1551 

-25.5799 



0.0000 

0.0000 



G = 




-0.4186 

-2.2548 

0.1121 

-1.2376 

-0.2107 

0.7058 

-0.0254 

1.6765 
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Table Bl. Continued 


a — 35° 

K - 

11.7391 

-58.3414 

1.8687 

-17.7790 

0.0000 

G = 

0.1880 

-0.3412 

a = 40° 

K 

13.9902 

-70.6014 

0.3247 

-19.3750 

0.0000 

G - 

0.1116 

-0.3338 

a — 45° 

K = 

10.6672 

-83.7570 

-6.8896 

-20.8578 

7.3404 

G = 

-0.4819 

0.1811 

a = 50° 

K = 

-11.8843 

-92.7368 

-11.4007 

-22.4435 

10.6763 


-33.3489 

-70.4668 

-14.4658 

-24.7922 

0.0000 


-2.2674 -0.2291 

0.7287 -1.3296 


-46.6544 

-56.1793 

-18.6618 

-22.5693 

0.0000 


-1.7216 0.0175 

0.3032 -2.8258 


-61.5933 

-38.0300 

-16.0345 

-20.3934 

-8.0219 


-0.8720 -0.6353 

0.3329 -0.2938 


-76.7698 

-22.1658 

-12.8837 

-18.4057 

-13.8889 


-0.9062 

1.3411 


-0.4289 

1.2669 


-1.0400 

1.5237 
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Table Bl. Concluded 


G = 


—0.1495 

-1.0156 

-0.2759 

-0.9231 

-0.0301 

0.4716 

0.2016 

1.3646 

a = 55° 




K = 




-29.1475 

-83.0683 



-89.7036 

-0.6993 



-2.1824 

-1.6032 



-25.2961 

-17.2914 



11.6376 

-18.5082 



G = 




-0.1753 

-0.8326 

-0.0765 

-0.7941 

-0.0716 

0.4357 

0.1557 

1.5059 

a = 60° 




K = 




-37.1976 

-83.9399 



-66.6277 

0.3374 



0.0000 

0.0000 



-27.7290 

-15.5958 



13.9318 

-27.0141 



G = 




0.4204 

-0.5479 

-0.3657 

-0.1400 

0.0113 

0.3407 

0.2077 

1.3010 
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